Exciton condensation and charge fractionalization in a topological insulator film 
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An odd number of gapless Dirac fermions is guaranteed to exist at a surface of a strong topological 
insulator. We show that in a thin-film geometry and under external bias, electron-hole pairs that 
reside in these surface states can condense to form a novel exotic quantum state which we propose 
to call 'topological exciton condensate' (TEC). This TEC is similar in general terms to the exciton 
condensate recently argued to exist in a biased graphene bilayer, but with different topological 
properties. It exhibits a host of unusual properties including a stable zero mode and a fractional 
charge ±e/2 carried by a singly quantized vortex in the TEC order parameter. 



Introduction. — Recent advances in studies of band in- 
sulators with strong spin-orbit coupling revealed the ex- 
istence of new topological invariants that characterize 
these materials [1 . Among the three-dimensional time- 
reversal (T) invariant insulators, the most interesting 
phase implied by this classification is the "strong" topo- 
logical insulator (STI), characterized by gapless fermionic 
states residing at its surface with an odd number of topo- 
logically protected nodes. These gapless states exhibit 
linear dispersion and behave as massless Dirac fermions 
familiar from the physics of graphene. Several real ma- 
terials have been identified as STIs in pioneering experi- 
ments [21 13] completed shortly after the theoretical pre- 
dictions f4j. These rapid developments give hope that 
the new state of quantum matter realized in STIs might 
be relatively common in nature and raise the prospects 
of future practical applications. 

The existence of an odd number of Dirac fermions leads 
to a number of exotic properties associated with surfaces 
of a STI. These include an exotic superconducting state 
induced by a proximity efi^ect that supports Majorana 
fermions [5] , a T-breaking phase which exhibits a frac- 
tional quantum Hall effect and an unusual 'axion' 
electromagnetic response [SJ [7] . 

The wealth of exotic phenomena listed above stems 
from the possibility of inducing various types of mass 
terms in the otherwise massless Dirac fermion states at 
the surface of a STI. In this Letter we introduce and 
study a new type of mass gap that can be induced by a 
Coulomb interaction between the surface states of a thin 
STI film and can be characterized as a 'topological' exci- 
ton condensate (TEC). The idea is motivated by recent 
proposals to realize an exciton condensate in a symmetri- 
cally biased graphene bilayer [1113] • We argue below that 
TEC in an STI film might be more easily realized than in 
graphene and is a different, genuinely topological phase, 
distinguished by the presence of a zero-energy mode and 
fractional charge associated with its vortices. 

Consider a film made of a STI placed inside a capacitor 
as in Fig. [l^. Imagine for simplicity that each surface 
harbors a single Dirac cone with the chemical potential 
/i initially tuned to the neutral point = 0. When the 




FIG. 1: (Color online) a) Schematic of the proposed device, 
b) The exciton condensate effectively joins the surfaces of the 
STI film resulting in toroidal topology. Arrows illustrate the 
magnetic field distribution of a planar monopole representing 
a vortex in the effective theory. 

capacitor is charged the Fermi levels in the two layers 
move in the opposite direction, creating a small electron 
Fermi surface in one layer and a small hole Fermi surface 
in the other. For arbitrarily weak repulsive interaction 
such a system will form an exciton condensate which may 
be pictured as a coherent liquid of electron-hole pairs 
residing in different layers. 

In what follows we use a simple model for the sur- 
face states to show how exciton condensation can be in- 
duced by the interlayer Coulomb interaction. By exam- 
ining this model we then deduce some interesting proper- 
ties of the underlying TEC. Specifically, we demonstrate 
that an isolated singly quantized vortex in the complex 
scalar order parameter characterizing TEC contains a 
zero mode and carries topologically protected exact frac- 
tional charge ±e/2. We put our findings in the context of 
axion electrodynamics, the low-energy effective theory of 
STIs, and discuss prospects for experimental realization 
and detection of the predicted phenomena. 

Model. — The gapless states associated with the two 
surfaces of the biased STI film can be described at low 
energies by a Dirac Hamiltonian [S] 

H = ^ il}l{vi(T --p- ^i)i)i + Unin2, (1) 
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where V'z = (qTi'^u)'^ denotes the fermion spinor in sur- 
face layer I, cr — (ax, cry) is the vector of Pauli matrices 
in the spin space, p = — iV, and vi = (— repre- 
sents the Fermi velocity, assumed to be opposite at the 
two surfaces and ni = ipjipi- We take in = fi + {—lYV, 
with /i the intrinsic chemical potential and V the external 
bias. The last term in ([T|) describes the short-range part 
of the interlayer Coulomb potential, a sufficient minimal 
interaction for the formation of EC [10] . Also, we assume 
the film to be sufficiently thick so that any direct hop- 
ping of low-energy electrons between the surfaces can be 
neglected. 

To describe the exciton condensation we decouple the 
interaction term in H using a matrix-valued order pa- 
rameter M = [/(V'lV'l) • The expectation value is taken 
with respect to the mean-field Hamiltonian 

Hmf = Ho + {i^\M^2 + h.c.) + iTr(AftM), (2) 

where Hq denotes the kinetic term in (jlj. At this point 
it is useful to organize the Fermi fields into a single 4- 
component spinor ^ = (■0i,'02)^- We can write Hmf = 
^'^H^' + ^Tr(M^M) with a 4 x 4 matrix Hamiltonian 

f v(T-p-V M \ 

Various forms of matrix M describe different possibil- 
ities for the TEC order parameter. When /i is close to 
zero an order parameter that opens up a gap in the ex- 
citation spectrum will be favored because it leads to an 
overall reduction in kinetic energy. In the uniform system 
this will occur only for M diagonal in the spin space, i.e. 
M = ml with m a complex constant, since the part of 
Ti. proportional to m then anticommutes with the kinetic 
term. The spectrum of Ti. then contains four branches 
and reads 

Ekas = -A* + ay^{v\k\ + sV")2 + |m|2, a, s = ±l. (4) 

Physically, this form of the matrix M implies non-zero 
expectation values (ci|C2|) = {ciic^^) — m/U. We also 
note that this is the only choice of the order parameter 
that leaves the Hamiltonian T-invariant. 

Exciton condensate. — When /i = and V =/= the ex- 
citon instability of our Hamiltonian ([T]) is formally equiv- 
alent to the Cooper instability in a metal and occurs at 
infinitesimal coupling U. To see what happens when fi is 
slightly detuned from zero, we consider the gap equation 
for TO (now assumed to be real) obtained by minimizing 
the ground state energy Eg = X^kas Ei^as + Nm'^ /U with 
respect to to. This reads 

TO 1 

kas 

where the prime denotes a sum over the occupied states 
{Euas < 0) only. For general values of fi, V and U the gap 



equation must be solved numerically. One can, however, 
extract the value Uc of the critical coupling beyond which 
TEC is formed. In the experimentally relevant regime 
/i ^ <C A, where A denotes the high-energy cutoff of 
the order of bandwidth, we find 

C/..4A(l + ^ln^) \ (6) 

The critical coupling is large, of the order of bandwidth, 
unless <C l^e""^/^, in which case Uc becomes small 
and eventually reaches zero when /i 0. In this limit 
the gap equation can be solved explicitly to obtain to w 
2^/VAe~-^ 1^^ . In order to achieve exciton condensation 
for a given coupling strength V it is essential to tune /i 
as close to zero as possible and apply high bias voltage 
V . Henceforth we consider only the = situation. 

Because of its exponential dependence on the coupling 
strength it is difficult to give a truly quantitative estimate 
of TO and the relevant TEC transition temperature 
for a realistic STI film. In the context of the graphene 
bilayer the estimates of range from sub-Kelvin up to 
the room temperature, depending on the approximation 
employed [8j [9l [12] . Although we do not attempt such a 
quantitative analysis here we note that the situation in 
STI might be quite similar. On the one hand the intrinsic 
energy scales in known STIs are somewhat smaller than 
in graphene. On the other hand, screening is known to 
reduce the mean- field by a factor ~ , where M 
denotes the number of surface Dirac modes. Formation 
of TEC in an STI film made from Bi2Se3 will be therefore 
aided by the fact that this material exhibits A/" = 1 [2] 
compared to A/" = 4 in graphene (due to valley and spin 
degeneracies) . 

Vortex zero modes. — In the following we adopt the 
point of view that, based on the above analysis, formation 
of TEC is likely to occur under experimentally achievable 
conditions and focus on its unique properties. To this end 
it is convenient to write the Hamiltonian (|3| in a more 
customary form using the Dirac matrices in the Weyl rep- 
resentation, 7j = iT2 (8) cTj, j = 1,2,3, 7o = Ti (g) 1 and 
75 = — i7o7i7273 = T3 (g) 1, where tj are Pauli matrices 
in the layer space. We obtain 

H = 70 {liPx + l2Py + V^7o75 + iTOle"'"*"^^) , (7) 

where we have set v = \ and used a polar representation 
TO — \m\e^^ of the complex TEC order parameter. 

The Hamiltonian (|7| formally coincides with the one 
used to describe the effect of EC near one of the val- 
leys of the biased graphene bilayer system [TU]. That 
work established existence of an exact zero mode of the 
Hamiltonian ([7| in the presence of a singly quantized 
vortex in the EC order parameter to — toqc*''', with 
(fi the polar angle. The zero-energy eigenstate has the 
form ^0 = {L9,i9\-irf with / = Ae-^-^ MVr), 
g = —iAe^^e "^'''^Ji{Vr) and A the normalization con- 
stant. In graphene, valleys always come in pairs. The 
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fermionic zero modes are thus doubled and split due to in- 
tervalley scattering. Because of this, no exact zero modes 
survive in the graphene bilayer. 

In a STI, by contrast, there is always an odd number 
of valleys associated with the surface [1] . When M — \, 
as in Bi2Se3, the zero mode attached to a singly quan- 
tized vortex will remain exact, as long as the vortex stays 
well separated from other vortices or system edges. This 
finding, along with the fractional charge discussed below, 
constitutes the key universal difference between the STI 
exciton condensate and other proposed exciton conden- 
sates and is the main result of this work. 

For general odd A/" > 1 we expect M —1 zero modes to 
split symmetrically around zero energy while the remain- 
ing zero mode will persist. This conclusion follows from 
the property 72^^*72 = 'H which together with 7! = — 1 
implies the spectral symmetry around the zero energy of 
eigenstates of TL. Specifically, for each eigenstate ^'^j of 
energy E there exists an eigenstate fi^ e with energy —E. 
Here O = 72 -ft^ is an antiunitary operator and K denotes 
complex conjugation. It also holds that {^,0.} = 0. 

Fractional charge. — A localized zero mode in a 
particle-hole symmetric system is known to carry a frac- 
tional charge ±e/2 [13l [HI [15l [16] . Thus, we expect our 
vortices to be fractionally charged. To see how this oc- 
curs in the present system and also to deduce some of its 
other interesting properties, let us consider an operator 
O, represented by a constant 4x4 Hermitian matrix, act- 
ing in the space of wavefunctions ^'£;(r). Following [TP 
consider now the quantity 



= Y.^^'E\0\^E)^\Ntr{0), 



(8) 



where N is the total number of quantum states in a suit- 
ably regularized system (e.g. on the lattice and in finite 
volume). The last equality in Q follows from the com- 
pleteness of states. We may write 



= 




O|^£) + (^'o|0|^'o) (9) 



The spectral symmetry generated by and its an- 
tiunitarity, expressed as {Q.^i\Q.^2) = (^'i|^'2)*, vca- 
Ply Ei=;>o(*s|0|vl/B) = S^<„(vI/^|(f]-iOf])t|vI/^). If O 
furthermore commutes with then the last term becomes 



simply ^E<o{^ e\0\^ e) and we can combine Eqs. (8][9) 
to obtain 



Y,{^e\0\^e) = 



E<0 



N 



tr(0) - (^'olOI^'o) 



(10) 



The expectation value of an observable represented by 
a constant 4x4 matrix that commutes with f2, taken 
over all occupied negative-energy eigenstates of Ti., is de- 
termined solely by the value of tr(0) and the zero-mode 
eigenstate of 7i. For an infinite system in continuum Eq. 
( 10 1 will be useful for quantities independent of N\ this 



occurs when O is traceless or else for quantities repre- 
sented as differences so that iVtr(O) drops out. In such 
cases the expectation value only depends on the zero 
mode and its value is expected to be robust. Specific 
examples follow below. 

The charge operator is represented by a 4 x 4 unit 
matrix Oq — el. The charge bound to a vortex can be 
expressed as 

Qv^eY, {{^E\t\^E)l - {^E\t\^E)o), (11) 

£;<o 

where subscripts 1 and refer to the state with one and 
zero vortices respectively. Using Eq. ( 10 1 we find Qv — 



(^'o|^'o)i — —§7 as expected. We note that Eq. (Ill 



assumed the zero mode to be unoccupied; if we occupy 
it by an electron then the vortex charge becomes +e/2. 

Other quantities of interest include the spin S and the 
axial charge carried by the vortex, defined as the 
charge difference between the layers. These are repre- 
sented by matrices Os = 570775 and Ogs — 675. Unfor- 
tunately these anticommute with U, making Eq. (10) in- 



applicable. A quantity that can be calculated is the inter- 
layer spin polarization AS, represented by Oas = \lol- 
A straightforward calculation shows that, in the pres- 
ence of a vortex, {/S.Sx) — {^Sy) — while (AS^) varies 
smoothly from ^ to as we tune V/m from to infinity. 
This implies that the vortex carries a fractional value of 
spin polarization between the layers. We note, however, 
that due to the entanglement of the spin and momenta 
in the STI, the total spin operators in the TEC or the 
STI surfaces are not sharp. 

The symmetry generated by f2 is a combination of T 
and spatial parity V. The latter will be broken in the 
presence of non-magnetic impurities and the zero mode 
will no longer be exact. However, the following general 
argument shows that the fractional charge remains pre- 
cisely quantized as long as the bulk is gapped and T is 
preserved. 

Axion electrodynamics. — As demonstrated in Refs. [5J 
[7] the response of STI to external electromagnetic field is 
that of an 'axion' medium [17! and can be mathematically 
implemented by adding a term 



AC, 



2TT he 



B E 



(12) 



to the usual Maxwell Lagrangian. An ordinary insula- 
tor has = while the STI exhibits 9 = ir, the two 
values permitted by the time-reversal symmetry. When 
a surface of a STI is gapped by a T-breaking perturba- 
tion, such as an applied magnetic field, 9 varies smoothly 
between tt and 0. As a result the surface behaves as a 
quantum Hall fluid [6l [17] . 

As noted above, TEC does not break T. Instead, the 
TEC order parameter effectively identifies the opposite 
surfaces of the film by allowing, at the mean-field level of 
the Hamiltonian ([2]), the electrons to hop between them. 
Dirac fermions at the surfaces are gapped without vio- 
lating T by a perturbation that effectively removes the 
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surfaces as illustrated in Fig. [TJd. Because the TEC order 
parameter is in general complex the electrons hopping be- 
tween the surfaces may acquire a nontrivial phase. This 
is easily included by postulating twisted periodic bound- 
ary conditions with a twist equal to the phase of the 
TEC order parameter x- The situation is easiest to vi- 
sualize in the context of a lattice model of STI where 
electrons traversing the bonds that connect the two sur- 
faces acquire a phase %. This extra phase can be also 
viewed as resulting from an electromagnetic vector po- 
tential 5 A = zx(<i>o/a) localized in the layer between the 
surfaces, with $0 — hc/e the flux quantum and a the 
lattice spacing. 

In a vortex configuration with xi^) = t^s it is easy to see 
that 5A(r) has the form of a planar monopole: the mag- 
netic field (5B radiates outward from the vortex center in 
the plane of the surface as illustrated in Fig. [T]d. The 
total fiux is (f>o- This vector potential is fictitious in the 
sense that (5B cannot be detected by an outside probe. 
To electrons, however, SA is indistinguishable from the 
real vector potential and must be included in Eq. ( 12 ) 



when evaluating the response of the STI. A unit mag- 
netic monopole in a '0- vacuum' described by Eq. ( 12 ) 
is known to carry electric charge — e(0/27r -I- n) with n 
integer [18]. Applied to TEC this gives vortex charge 
— e( 1/2-1- n), consistent with our finding of the fractional 
charge ±e/2 bound to the vortex. 

Outlook and open questions. — Recent advances in ma- 
terials engineering and fabrication give hope that the ex- 
otic state of matter identified in this work can be achieved 
and probed in the near future. From an experimental 
point of view, there are several significant advantages 
of the proposed phase. First, we believe that making 
an exciton condensate between the surfaces of a single 
film, rather than from two 2D materials (e.g., graphene) 
with an insulator between them, is easier because it does 
not require creating a pinhole-free insulating layer with 
defect-free junctions to the two 2D materials. Second, at- 



taching leads to a surface of a film should be considerably 
easier than attaching leads to graphene. Once the leads 
are in place it should be straightforward to identify the 
onset of the exciton condensation in a transport measure- 
ment as demonstrated in recent literature [191 . Third, we 
note that the existence of the zero-energy Dirac fermion 
can be detected in the same way as a standard midgap 
impurity state by optical methods or in a careful trans- 
port measurement. In the TEC phase a transport mea- 
surement will reveal a gap at low temperature. As the 
temperature nears the transition temperature, conduc- 
tion becomes dominated by the charges bound to vortices 
and is proportional to the number of thermally excited 
vortices. Fractional charge can be probed, at least in 
principle, by the shot noise analysis of resistivity [20 . 

From a theoretical point of view TEC in STI film is 
fundamentally interesting for several reasons. To our 
knowledge, TEC is the first example of a new symmetry- 
breaking phase enabled by the special properties of topo- 
logical insulators. It differs from the superconducting 
state generated at the surface by proximity effect [5] be- 
cause in that case there is no new symmetry breaking. 
It differs from the ordinary exciton condensate, which 
is in the same universality class as a ''He film, because 
the zero mode attached to a vortex is stable. In the 
case of graphene bilayer the zero modes are not protected 
due to inter-valley scattering. In this respect as in sev- 
eral others, topological insulators allow the realization of 
physics that is spoiled in graphene by inter-valley scat- 
tering. The consequence is a distinct low-temperature 
phase of matter with fractionally charged topological ex- 
citations whose exchange statistics presents an interest- 
ing open question. 
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